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Abstract 

Recently, Zhang et al. [PRA, 75, 062102 (2007)] extended Kieu's interesting work on the quan- 
tum Otto engine [PRL, 93, 140403 (2004)] by considering as working substance a bipartite quantum 
system AB composed of subsystems A and B. In this paper, we express the net work done Wab by 
such an engine explicitly in terms of the macroscopic bath temperatures and information theoretic 
quantities associated with the microscopic quantum states of the working substance. This allows 
us to gain insights into the dependence of positive Wab ° n the quantum properties of the states. 
We illustrate with a two-qubit XY chain as the working substance. Inspired by the expression, 
we propose a plausible formula for the work derivable from the subsystems. We show that there 
is a critical entanglement beyond which it is impossible to draw positive work locally from the 
individual subsystems while Wab is positive. This could be another interesting manifestation of 
quantum nonlocality. 

PACS numbers: 03.65.Ud, 07.20.Pe 
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Heat engines are devices that extract energy from its environment in the form of heat 
and do useful work. At the heart of every heat engine is a working substance, such as a 
gas-air mixture in an automobile engine. The operation of the heat engine is achieved by 
subjecting the working substance of the engine to a sequence of thermodynamic processes 
that form a cycle, returning it to any arbitrarily selected state. Quantum heat engines, in 
contrast, operate by passing quantum matter through a closed series of quantum thermo- 
dynamic processes [l|]. For instance, Kieu 

as 

introduced a class of quantum heat engines 
which consists of two-energy-eigenstate systems (qubits) undergoing, respectively quantum 
adiabatic processes and energy exchanges with heat baths at different stages of a cycle - 
the quantum version of the Otto cycle. Recently, Zhang et al. (4] extended Kieu's work by 
considering the quantum Otto engine with a two-qubit (isotropic) Heisenberg XXX chain 
as the working substance. The chain is subject to a constant external magnetic field. The 
purpose of their paper is to analyze the effect of quantum entanglement on the efficiency of 
the quantum Otto engine. This is an important and intriguing development since it brings 
together concepts from quantum mechanics and thermodynamics - two seemingly different 
fundamental areas of physics. 

Entanglement is a property associated with the state of a composite quantum system 
made up of at least two subsystems. It is a nonlocal correlation between quantum systems 
that does not exist classically. Therefore, it is imperative to raise the following questions. 
What is the relationship between the net positive work derivable from the subsystems and 
that from the total system? What is the role of entanglement in this relationship? In 
order to obtain plausible answers to these questions, we need a means to calculate the work 
derivable from a subsystem. In this paper, we propose information theoretic answers |||. 
As a concrete example, we consider the two-qubit XY model [6j. The Hamiltonian for the 
two-qubit XY chain in an external magnetic field B m along the z axis is given by 

H = -(1 + i)J<j\ ® a\ + ^(1 - 7) Ja\ ®o 2 B + -B m (a\ <g> a° B + a° A ® 0%), (1) 

where er° is the identity matrix and a l a (i = 1, 2, 3) are the Pauli matrices at site a = A, B. 
The parameter — 1 < 7 < 1 measures the anisotropy of the system. (I + 7) J and (1 — 7) J are 
real coupling constants for the spin interaction. The chain is said to be antiferromagnetic 
for J > and ferromagnetic for J < 0. Here, we consider J > 0. 

To set the stage, we describe the four quantum thermodynamic processes of the quantum 
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Otto cycle. In the following, we consider as working substance a bipartite quantum system 
consisting of two subsystems, A and B, with Hamiltonian H = J2i Ei\^ 1 ) ab(^ 1 \- Here, Ei 
and \^ % ) ab are respectively the eigenvalues and eigenvectors of H. For the two-qubit XY 
model, we have 

l^ 1 )^ = 1 [(B + B m )\00) A B + lJ\n)AB}, 

'(B + B m ) 2 + 1 2 J 2 



\^ 2 )ab = ^=[\01)ab + \10) ab], 
\^)ab = -^[|01)ab-|10)ab], 

\^) ab = i 1 =[{B - B m )\m) AB - iJ\±1)ab\, (2) 

>(B - B m f + 7 V 2 



with E\ = — E4 = B and Ei = —E3 = J, where B = yB^ + ^J 2 . Furthermore, we 
assume that the system is allowed to thermalize with the heat baths in processes 2 and 4. 
Specifically, process 4 brings the system to its initial quantum state given by the density 
operator 

P% = Y,P^ %1 )ab{V% (3) 

i 

with pn = exp(—En/kTi)/Zi and Z\ = J2i ex P( — En/kTi). That is, the system is initially 
in thermal equilibrium with a heat bath at temperature T\. For the two-qubit XY model, 
we have E iX = E { and \^ i1 )ab = \^ % ) ab with J = J x . 

1. The system is isolated from the heat bath and undergoes a quantum adiabatic process, 
with for instance J changing from J x to J 2 - Provided the rate of change is sufficiently 
slow, Pii's are maintained throughout according to the quantum adiabatic theorem [7]. 
At the end of process 1, the system has the probability p a in the eigenstate \^ i2 )ab- 
According to Kieu 

as 

, an amount of work is performed by the system, but no heat 
is transferred during this process. 

2. The system is brought into some kind of contact with a heat bath at temperature 
T2 < T\. After the irreversible thermalization process, the quantum state of the 
system is given by the density operator 

pS = EpoI*°>^<*°I. (4) 
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where pi 2 = exp(— E i2 /kT 2 )/Z 2 and Z 2 = X]; ex P( — Ei 2 /kT 2 ). Here, for the two-qubit 
XY model, we have E& = Ei and \^ i2 )ab = ab with J = J 2 - It follows from 

as 

that only heat is transferred in this process to yield a change in the occupation 
probabilities, and the heat transferred is given by 

Q 2 = Y^E i2 (pi2 -Pn)- (5) 

i 

3. The system is removed from the heat bath and undergoes a quantum adiabatic process, 
with for instance J changing from J 2 to J\. At the end of process 3, the system has 
the probability pi 2 in the eigenstate \^ i1 )ab- An amount of work is performed on the 
system, but no heat is transferred during process 3. 

4. The system is brought into some kind of contact with a heat bath at temperature 
T\. After the irreversible thermalization process, the quantum state of the system is 
returned to the initial one in Eq.Q. The heat transferred in process 4 is given by 

Q 4 = J2 E a{pa -Pa)- (6) 

i 

From the first law of thermodynamics, the net work done by the quantum Otto engine 
is Q,H 

W AB = Q2 + Q4 

= ^2(E a -Ei2)(pa-Pa). (7) 

i 

Substituting = —kTj\og(pijZj) into Eq.fljTJ) and setting the Boltzmann constant k = 
log 2 e, we obtain 

W AB = (T, - T 2 ){S[p%) - S[pf B }} - TtflpaWpa] - T 2 H\p a \\p i2 }. (8) 

Here, SVab] is the von Neumann entropy of the quantum state p^ B , and if[p i2 |bii] and 
H\pn \\Pi2} are the relative entropies of p i2 to pn and pn to p i2 respectively It follows from 
the non-negativity of the relative entropy that to derive positive work we not only have to 
demand T% > T 2 but also 5"[p^] > S[pa B ] such that if[pj 2 |bii] an d H\pn\ \p i2 ] are not too 
large. This is true regardless of whether AB is a single quantum system or one composed 
of two or more subsystems. We shall illustrate this after the following remark. 
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We note that iJ[p i2 |bii] = ^[PabII/mbL the quantum relative entropy of p^ B to p^ B 
if they share the same eigenstates. Similarly, we have if^l |p i2 ] = £[Pab| \Pab\- This is 
the case for the two-qubit XY model when we let B m = rjj, with rj some fixed constant. 
We assume this to hold from here on without loss of generality. It follows immediately 
from Eq.([7|) that Wab is directly proportional to (Ji — J 2 ) for any given 7 and 77. Another 
consequence is that p^'s and p^'s depend solely on J\/T\ and J 2 /T 2 respectively. It then 
follows from Eq.® that Wab = when J 2 /T 2 = Ji/Ti. So, in order to have positive Wab, 
we require (T 2 /Ti) J x = J m m < J 2 < Ji- 
lt is clear, from Eq.([8]), that the condition J min < J 2 is one on the quantum states p^ B 
(2) 

and Pab- ^ * s thus natural to express this condition in terms of quantities that describe 
the two-qubit states. Here, we recall the quantum mutual information between the two 
subsystems A and B, 2^(A : B) = %$' } ] + S\pf]- S\p%] fl. This is usually used 
to measure the total correlations between A and B [9]. Hence, in order to derive positive 
Wab we demand that 1^ 2 \A : B) > X <yl \A : B) = X^ n . If the states are not separable, 
we may require that C[p^b] > £[/°ab] = QxLj where C[p^] is the Wootters concurrence 

js] associated with Pab- This measures the amount of entanglement or nonlocal quantum 

(2) (2) 

correlation between A and B. Given Ti and J\, there is therefore a minimum or C^ in 

below which the quantum Otto engine does not yield any positive work. 

Equation ([8]) gives explicitly the dependence of Wab on the bath temperatures Ti and 

T 2 , and on the quantum states p^ B and p^ B . Consider some T 2 , J min < J 2 < Ji and kT 2 , 

(2) 

«;J 2 , which yield identical p\ B . Here, k is a positive constant such that kJ 2 > J\. Therefore, 
Wab is positive in the first case but negative in the latter since nJ 2 > J\- This difference 
is obviously due to the bath temperatures. Our focus here is on analyzing the dependence 
of positive Wab on the quantum properties of the states. Therefore, for given Ti and J\, it 
suffices to chose an appropriately small T 2 such that as J 2 is increased from J m \ n to J\, we 
have states Pa B , with quantum mutual information going from Z^ n to sufficiently close to 

/OA 

the maximum possible of 2 (or with concurrence going from C^J^ to sufficiently close to the 
maximum possible of 1). 

Equation ([HD also provides quantum information theoretic insights into the condition 
J 2 < Ji. Suppose 7 = 0.4, 77 = 0.3, Ti = 1000, J x = 8, and T 2 = 0.1. Then, p n « p 21 « 
P31 ~ Pn ~ 0.25 and p^' B approximates the density operator of a maximally mixed state with 
neither classical nor quantum correlations. We also note that T^ 2 \A : B) m 2 and C[p^g] ~ 1 
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FIG. 1: Net work done Wab by the quantum Otto engine with a two-qubit XY chain as the 
working substance plotted vs J 2 , for 7 = 0.4, r\ = 0.3, T\ = 1000, J\ = 8, and T2 = 0.1. 
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FIG. 2: The figure on the left hand side shows the behaviour of (Ti - T 2 ){S[p { } B ] - S[p%]} (gray) 
and TiH\pi2\\pii\ (black), while that on the right hand side shows the behaviour of T 2 H\pn\\pi 2 }. 

when J 2 = 8, satisfying the above sufficient condition. Figure 1 shows the dependence of 
Wab on J 2 . Wab increases from zero at J 2 = J m - m = 8 x 1CT 4 to a maximum W max 10.3695 
at J 2 = J max ~ 0.575065, after which it decreases from W max to zero at J 2 = J±. Increase in 
J 2 yields p^ B with (T 1 — T 2 ){S[p^ B \ — S\p^ B \\ and TiH\p i2 \\pii] approaching the approximate 
maxima 1999.81 and 1988.48 respectively, but with T 2 H[pn\ \pi 2 ] that increases monotonically 
with J 2 (see Fig. 2). At J 2 = 8, T 2 H[pii\\p i2 ] exactly equals the difference between the 
maxima. And, Wab becomes negative if J 2 is increased beyond J 2 = J\. Intuitively, one 
would expect to draw more work as the von Neumann entropy of p^ B differs more from 
that of p^ab'i for instance, when p^ B becomes increasingly correlated. This is indeed the case 
before this difference results in significant increases in the relative entropies enough to cause 
Wab to decrease to zero. For J min < J 2 < J max , Wab increases with increasing quantum 
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FIG. 3: Quantum mutual information in the two-qubit XY chain plotted vs J2, for 7 = 0.4, 

T) = 0.3, and T2 = 0.1. At J2 = 0, it is zero since the qubits do not interact. It then increases with 
increasing J2 to a maximum of 2 log 2 where the correlation is completely quantum in nature. 
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FIG. 4: The state p^ B is separable when J m - m < J2 < 9.31358 x 10~ 2 . Otherwise, it is entangled - 
with the concurrence approaching one as J2 increases since P32 - > 1 and p\2 ~ P22 ~ P42 — ► with 
increasing J2. 

mutual information (compare with Fig. 3). During this increase the correlation changes 
from purely classical to a mixture of classical and increasingly quantum ones (see Fig. 4). 
Beyond J 2 = J max when C{p%] « 0.903444, W AB be gins to decrease. It is not obvious at 
this stage what precisely is the role of entanglement. We shall further elaborate on this after 
we propose the following definition for the work derivable from the individual subsystem. 

The subsystems A and B are clearly being subjected to exactly the same quantum Otto 
cycle that the composite system AB has undergone. The quantum states of A and B that 
correspond to Eqs.flB]) and (j4j) are also well-defined. Namely, = tr Bp ab an d Pb = 
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FIG. 5: Net work done wa by the quantum Otto engine with one qubit of a two-qubit XY chain 
as the working substance plotted vs J2, for 7 = 0.4, 77 = 0.3, T\ = 1000, J\ = 8, and T2 = 0.1. 

ti aPab- Hence, inspired by Eq.flBD, we define 

w a = {T x - T 2 ){S[p«] - S[p^}} - Tx^ll^] - T t H[qg>\\qg>], (9) 

where a = A or B, q^fs and s are the eigenvalues of p^ and p B respectively. Now, 
consider p^ B approximately maximally mixed and p^ B sufficiently close to a maximally 
entangled state, then both an d p£) will be almost identical to the maximally mixed 
state. According to Eq.(jHJ), the work derivable from a subsystem undergoing the quantum 
Otto cycle in this case will be extremely close to zero. For the two-qubit XY model, p^ and 
Pb are identical. Figure 5 shows wa as a function of J2. Wa increases from zero at J2 = J m in 
to a maximum w max « 0.231128 at J 2 = j max ~ 0.166289, when C\p^ B ] « 0.316188. After 
T2 = jmax, u^a decreases from u> max to zero at J 2 = jcnt ~ 1.24252 when C[p^] = C^ t ~ 
0.9964. Therefore, Eq.Q quantitatively describes the above consideration. We propose 

here that wa {wb) is the work derivable from the subsystem A (B). Therefore, given T\ and 

(2) 

Ji, there is a critical concurrence Q rit above which no positive work can be drawn locally 
from each subsystem. Increasing J 2 beyond j crit , wa becomes negative. This is because pj± 
becomes more random than p^ . 

In conclusion, we have expressed, in Eq.flS]), the net work done Wab by a quantum Otto 
engine explicitly in terms of the macroscopic bath temperatures (Ti and T 2 ) and information 
theoretic quantities associated with the microscopic quantum states (pab an d Pab) °f the 
working substance - a bipartite quantum system. Equation (|Sj) inspires our proposal of a 
plausible formula, Eq.flJj]), for the work derivable from the subsystems A and B. For the 
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two-qubit XY chain, we show that p^ B must have quantum mutual information (A : B) 
or concurrence greater than that of a given p$ B to yield positive Wab- Equation (JHJ) 

also provides information theoretic insights into how Wab increases and then decreases with 
increasing 1^(A : B) or C[p%}. We show, using Eqs.flE]) and ([9]), that there exists a critical 
concurrence above which no positive work can be derived locally from each subsystem while 
Wab > 0. This could be another interesting manifestation of quantum nonlocality. 
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